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Abstract. Using LDA+GTB (local density approximation+generalizcd tight- 
binding) hybrid scheme we investigate the band structure of the electron-doped high- 
ly material Sni2- x Cc :Z ;Cu04. Parameters of the minimal tight-binding model for this 
system (the so-called 3-band Emery model) were obtained within the NMTO (N- 
th order Muffin-Tin orbital) method. Doping evolution of the dispersion and Fermi 
surface in the presence of electronic correlations was investigated in two regimes of 
magnetic order: short-range (spin- liquid) and long-range (antiferromagnetic metal). 
Each regime is characterized by the specific topologies of the Fermi surfaces and we 
discuss their relation to recent experimental data. 
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1. Introduction 



One of the most important questions in condensed matter is how the strong interaction 
between quasiparticles modify their properties and influence observable quantities. Non- 
Fermi-liquid behavior was found in many different substances, but a class of high-T c 
copper oxides attracts special attention during the last few decades. The unconventional, 
non-s-wave, superconductivity has a lot to do with it. While other players came to 
stage, like lamellar sodium cobalt oxides and Iron-based pnictide superconductors, only 
high-T c cuprates combine both strong electronic correlations and high values of critical 
temperature. 
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A key issue in a theory of high-T c superconductivity is the proper description 
of the low-energy electronic structure. Recent experimental results, mainly of angle- 
resolved photoemission spectroscopy (ARPES) [H El E] and measurements of quantum 
oscillations [U El E], provide a pattern to test various theoretical models and schemes. 
One of the approaches, proposed by some of the present authors, is the LDA+GTB 
hybrid scheme [7]. It was shown that the mean-filed theory within this scheme captures 
the most essential features of the doping-dependent evolution of the quasiparticle band 
structure and the Fermi surface [El [9]. 

A lot of theoretical and experimental efforts were concentrated on the hole doped 
compounds. Systems with electron doping, Re2-zCea;Cu04 (Re=Nd, Pr, Sm), present a 
counterpart to hole doped ones and a test for electron-hole asymmetry in Mott-Hubbard 
insulators. Recent ARPES data on the optimally doped, x = 0.14, Sm-based compound 
provide detailed information on the Fermi surface and the band dispersion in the vicinity 
of the Fermi level [10]. Similar study was reported for Nd-based compound by Schmitt 
et al. [H]. These results were confirmed independently by the measurement of the 
quantum oscillations [12]. Also, Park et al. [10] presented an explanation of the observed 
data based on the a/2 x s/2 spin-density wave (SDW) model. On the other hand, the 
high-energy electronic structure is found to be inconsistent with the SDW scenario. 
Moreover, the SDW model implies the weak or moderate electronic correlations and 
a Fermi liquid background, which is obviously not the case for the underdoped and 
optimally doped cuprates. Thus it is not a satisfactory scenario and a strong correlation 
effects should be taken into account. 

Here we present the investigation of electronic structure for the electron-doped 
high-T c material Sm2_ x Ce x Cu04 by means of LDA+GTB hybrid scheme. Parameters 
of a minimal generic tight-binding model for these systems (the so-called Emery model) 
were obtained within the NMTO method. Doping evolution of the band structure and 
the Fermi surface within this model in presence of strong electronic correlations and 
magnetic fluctuations were studied in the framework of the GTB method. 

The LDA+GTB electronic structure strongly depends on the underlying magnetic 
order. Though there is no Neel temperature in the optimally doped n-type cuprates, 
the antiferromagnetic (AFM) correlation length is extremely large (A ~ 400a) up to 
x = 0.17 [131 El [El HE] . Such correlation length makes magnetic behavior to be rather 
close to the long-range ordered AFM. That is why our spin-liquid description is unable 
to capture some details of the observed Fermi surface. On the other hand, there is 
a good agreement with the recent ARPES data once we assume the presence of the 
long-range order. 

2. Noninteracting band structure 

Here we will describe the noninteracting band structure and in the next Section introduce 
the electronic correlations within the LDA+GTB method. 

Sm 2 Cu04 system has body-centered tetragonal crystal structure with the space 
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Figure 1. Electronic structure of Sni2Cu04 obtained within LDA. Left panel presents 
total and partial densities of states; right panel shows band dispersions: thick (cyan) 
curves denote LMTO bands, while thin (black) curves denote NMTO bands. Zero 
corresponds to the Fermi level. 

group IA/mmm. Values of lattice parameters are a = 3.917Aand c = 11.899A. The 
atomic positions for different atoms are: Cu (0,0,0), Sm (0,0,0.35184) and two types 
of oxygens 01 (0,0.5,0), 02 (0,0.5,0.25) [IT]. Physically important Cu0 2 layers are 
constructed with 01 type oxygens. No apical oxygen is presented in this structure. 

We perform density functional theory band structure calculations within linear 
muffin-tin orbital basis set employing atomic sphere approximation in the framework 
of program package TB-LMTO-ASA v47 pi [191120]. In Fig. [1] results of our LMTO 
computations are presented. Left panel shows the total and the partial densities of 
states. Cu-3d and 01-2p states cross the Fermi level. In the right panel of Fig. [1] LMTO 
band dispersions are presented (thick curves). Note the Fermi level is crossed by just one 
antibonding hybrid Cu-3<i — 01-2p band of x 2 —y 2 symmetry. It is in agreement with the 
generic minimal tight-binding model for high-T c cuprates [21} [22]. Orbital basis for this 
model consists of C\x-"id x 2_ y 2 orbital and in-plane p x and p y oxygen orbitals. To compute 
corresponding model parameters iV-th order Muffin-Tin orbital method (NMTO) [23] 
was used. Necessary for NMTO expansion energies are schematically shown on the right 
side of Fig. [1] Obtained hopping parameters are listed in the Table d} and the single 
electron energies are E x 2_ y 2 = —2.322 eV and E Px = —3.708 eV. 

We assume that the values of Coulomb repulsion U and Hund's exchange Jh for 
Cu ions are doping independent and equal to 10 eV and 1 eV, respectively (see Ref. [7] 
for details). 
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Table 1. Parameters for the 3-band model obtained within NMTO method. Here x 2 
and p x ^ y denote the Cu-3d x 2_ y 2 and 0-p x ,y orbital indexes. All values are in eV. 
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3. LDA+GTB scheme 

Within the LDA+GTB method [7J the results of ah initio band structure calculations 
presented in the previous Section are used to construct the Wannier functions and 
to obtain the parameters of the multiband Hubbard-type model. For this multiband 
model the electronic structure in the strong correlation regime is calculated within the 
generalized tight-binding (GTB) method [2U [25j EE]. The latter combines the exact 
diagonalization of the model Hamiltonian for a small cluster (unit cell) with perturbative 
treatment of the intercluster hopping and interactions. After this step we end up with 
the GTB Hamiltonian. Depending on the order of perturbative treatment, we can 
formulate different approximations. 

As was shown before [7J, for undoped and weakly doped La 2 _ a; Sr : j;Cu04 and 
Nd 2 _ a; Ce :r Cu04 this scheme in the lowest order in hopping (Hubbard-I approximation) 
results in a charge transfer insulator with a correct value of the gap E ct and the dispersion 
of bands in agreement with the experimental ARPES data. 

We map the GTB Hamiltonian onto the effective t — t' — t" — J* model, where "*" 
denotes the three-site correlated hoppings, and study two regimes of AFM correlations. 
Namely, the spin-liquid phase with short-range AFM fluctuations and the long-range 
AFM metallic phase. We will describe AFM metallic phase using the Hubbard-I 
approximation that was shown to be in qualitative agreement with the Quantum Monte- 
Carlo results [27]. To study the spin-liquid phase, we use the same procedure as in 
Ref. [8] and go beyond the Hubbard-I approximation: (i) We solve the Dyson equation 
in the paramagnetic phase by means of the diagram technique for the Hubbard X- 
operators [28] and (ii) Obtain the coupled equations for the self-energy E(k, to), the 
strength operator P(k, u), and the spin-spin and kinematic correlation functions, (iii) 
We solve the coupled equations self-consistently and obtain a doping dependent Fermi 
surface and band structure. 

The t — t' — t" — J* model Hamiltonian is given by 

H t -j* = Ht-j + #(3), (1) 
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n) (n'\ are the Hubbard X-operators [29] acting on the Hilbert space of 
local states \n) = {0, a, —a = a}, Jf g = 2t"j g /U is the AFM exchange between two sites 
/ and g, U = E ct is the effective Hubbard repulsion determined by the charge transfer 
energy E ct « 2eV, t^ = 2t(cosk x + cosk y ) + At' cos k x cos k y + 2t"(cos2k x + cos2/c y ) is 
the Fourier transform of the hopping tf g , and t^ = 2t(cos k x + cos fey) + At' cos /c x cos & y + 
2t"( cos 2/c x + cos2k y ) is the Fourier transform of the interband hopping parameter tf g , 
Sf is the spin operator, e is the one-hole local energy, and /i is the chemical potential. 
The Green function in terms of the Hubbard X-operators is 

ctO\\ (l+x)/2 



GJk,u) 



X 



Oct 



X 



/U 2~tk — 



P 
Ik 



S(k) 



(2) 



Within our approximations [5], the strength operator P(k,E) is replaced by the 
occupation factor (1 + x)/2 and the self-energy S(k, ui) is frequency independent but 
preserve momentum dependence, 
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The spin-spin Cq and kinematic Xq correlation functions play significant role 
representing the short-range AFM fluctuations and the kinetic energy reduced by the 
correlation effects, respectively: 
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Energy spectrum is determined by the poles of the Green function (j2J) and Fermi 
surface is determined by the equation Sq — ji+ ^^t^ + ^—f-jf + S(k) = 0. 



4. Results and discussion 



The procedure of mapping the GTB Hamiltonian onto the effective model was described 
in detail in Ref. [7]. Following the same steps and using the parameters listed in Table [TJ 
we obtain the t—t'—t" — J* model with the following hoppings and exchange interactions: 
t = -0.59 eV, t' = -0.08*, t" = 0.15t, J = 0.92|t|, J' = 0.01|t|, J" = 0.02|t|, t = -0.74 
eV, t' = — 0.1 It, t" = 0.16t. Here, t, t', and t" are the interband hoppings through 
the charge-transfer gap, which determine the three-site hoppings and the exchange 
parameter, J = 2t 2 /E ct . Note that although the value of the nearest-neighbor exchange 
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Figure 2. Spin-liquid phase: Band structure and DOS (a) and the Fermi surface (b) 
for Sni2- 2 ;Cc a ;Cu04 within LDA+GTB method for different doping concentrations x, 
as indicated. In (a) zero corresponds to the Fermi level. 



J is quite large, the spin gap in the AFM phase will be determined not by this value 
alone, also there will be a contribution from the three-site hoppings. This contribution 
reduce the value of the spin gap as will be discussed later. 

In Fig. [2] we present results for t — t! — t" — J* model in the spin-liquid phase. 
At low doping, x = 0.03, due to the scattering on the magnetic fluctuations, the band 
structure possess local AFM symmetry in the vicinity of the (±7r/2, ±7r/2) points [see 
Fig. E£a)] and the Fermi surface has a form of four electron pockets around the (0, ±7r) 
and (±7r, 0) points. Values of the spin-spin correlation functions Cq are large enough for 
the similar topology to survive until x ~ 0.22, where a quantum phase transition with 
change of the Fermi surface topology takes place. After the transition, the Fermi surface 
at x = 0.22 has a form of a large hole pocket around the (n, ir) point and a small hole 
pocket around the (0, 0) point, which decrease in size with further electron doping. At 
x = 0.25 only one large hole pocket around the (tt, it) point is left. The quantum phase 
transition with the change of the Fermi surface topology was found experimentally in 
Nd-based compound [12], though at a different critical concentration. 

Note that the standard formulation of the Luttinger theorem does not work for 
the Hubbard fermions since the spectral weight of such fermion is determined by the 
strength operator, P(k, E) = F 0cr , and each quantum state contains 2F 0a = 1 — x 
electrons. A generalized Luttinger theorem for the strongly correlated systems [32] 
takes into account the spectral weight of each \k) state and the Fermi surface in our 
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Figure 3. Fermi surface "mapping" (set of equipotential cuts from the Fermi level to 
-0.3 eV below it) in the spin-liquid phase for x = 0.14 together with the reproduction 
of experimental ARPES data from Ref. [TU] in the lower right corner. Colorbar on 
the right shows the correspondence between shades of gray and energy from the Fermi 
level (in eV). 

Fig. [2] satisfies its completely. 

A comparison of the calculated Fermi surface in the spin-liquid phase and the 
experimental ARPES data [10] is shown in Fig. [3j Note the difference in the methods 
to obtain the Fermi surface "mapping" . We draw a set of constant energy cuts from the 
Fermi level down to -0.3 eV below it, while the experimental Fermi surface mapping is 
an integration of ARPES intensities over 30 meV energy window. The ARPES Fermi 
surface consists of three parts: two pockets around (0, it) and (ir, 0) points, and one 
elongated pocket around (tt/2,tt/2) point. One can immediately notice from Fig. [2(b) 
that in our spin-liquid theory the pocket around (7r/2,7r/2) point is missing; it does 
not appear even if one collects intensities from below the Fermi level, as seen in Fig. [3j 
Moreover, there are no features in the band dispersion, which could produce such pocket. 
Thus we conclude that our theory for the spin-liquid phase does not reproduce all details 
of the experimental Fermi surface. 

Since optimally doped Sm2_ x Ce x Cu04 is in the vicinity of the ordered AFM phase 
and the correlation length is extremely large (about 400 lattice constants) [IE] , we now 
investigate the band structure in the t — t' — t" — J* model assuming the long-range 
AFM order. The procedure is similar to Refs. [301 EI] i where the energy spectrum of 
the t — t' — t" — J model was obtained within the Hubbard-I approximation, but here 
we also take the three-site hopping terms into account. 

In Fig. H](b) we present results for the Fermi surface in the AFM phase of the 
t — t' — t" — J* model at x = 0.14 together with the experimental ARPES data. 
Evidently, there is a rather good agreement between both. We would like to mention 
that for lower concentrations our calculations result in decrease of pockets around the 
(tt/2, 7r/2) points and increase of pockets around the (n, 0) and (0, ir) points. Note that 
in the t — J model the spin gap is determined solely by the AFM exchange J [30]. Here, 
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Figure 4. AFM phase: Band structure (a) and the Fermi surface (b) for x = 0.14 
within LDA+GTB method. In the lower right corner of (b) we show the reproduction 
of the ARPES Fermi surface map from Ref. [10] . In (a) zero corresponds to the Fermi 
level. 



momentum dependence of the spin gap is proportional to t' cos k x cos k y (see Eq. (12) of 
Ref. [30]). The reason is that in the absence of spin fluctuations the hoping of a particle 
without spin flip processes possible only within the same spin sublattice. Because of the 
cos k x cos k y functional form the spin gap is maximal at (tt/2, 7r/2) point and minimal 
at (vr, 0) point as seen in Fig. H](a). Since the three-site hopping terms involving sites 
/, m, and g are proportional to tf m t mg / E ct , they also contribute to the spin gap; but, 
apparently, they decrease the gap value around (%, 0) point making it more anisotropic. 

Since we are making a mean-filed theory (though in a strong interaction limit) we 
can not address the question of the intensity distribution over the Fermi surface. This 
question was addressed earlier by different groups j33j G2]. Remarkably, their results 
on the Fermi surface contours for x ~ 0.14 are very similar to our's in Fig. H] in spite 
of rather different calculation schemes. This again emphasizes the fact that the AFM 
correlations are very strong in the optimally electron doped cuprates and they determine 
the quasiparticle dispersion and the Fermi surface. 

5. Conclusion 

We have shown that the experimentally observed Fermi surface topology can be 
explained within the LDA+GTB calculations for the long-range AFM spin background. 
On the other hand, our theory for the spin-liquid phase demonstrates only partial 
agreement with the ARPES Fermi surface due to the underestimation of the impact of 
magnetic scattering on the electronic structure. We conclude that the spin fluctuations 
are very strong in Smi.seCeo.uCuC^ and are closer to the long-range AFM fluctuations 
rather than to the fluctuations in the spin-liquid phase. Similar conclusion was drawn 
recently from the analysis of quantum oscillations in Nd-based electron doped cuprates 
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We would like to emphasize the significant difference between our picture for AFM 
order and one by Park et al. |10j . Park et al. provide a simple calculation based on 
the conventional SDW order i.e. the one based on a weak coupling approximation for 
the interaction. In the absence of the long-range order the ground state is metallic even 
at zero doping, x = 0. On the other hand, our approach allows to study the limit of 
large interaction and provides an insulating ground state at zero doping. This is essential 
difference since the underdoped cuprates belong to a class of strongly interacting systems 
and exhibit a Mott transition at a half filling, x = 0. More precisely, because of the 
copper-oxygen hybridization the cuprates shows the charge-transfer gap E ct at x — 0, 
but on the level of a single-band Hubbard model one can speak about a Mott-Hubbard 
effective gap U = E ct . 
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